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The use of simulated and controlled acoustic emission signals forms the basis of a 
powerful tool for the detailed study of various deformation and wave interaction processes in 
materials. We report the results of experiments and signal analyses of AE resulting from point 
sources such as various types of indentation-produced cracks in brittle materials and the growth 
of fatigue cracks in 7075-T6 aluminum panels. Furthermore, we describe our recent work 
dealing with the modeling and subsequent signal processing of an extended source of emission in 
a material. Results of the forward problem and the inverse problem are presented with the 
example of a source distributed through the interior of a specimen shown. 


INTRODUCTION 

The use of acoustic emission (AE) signals for real-time structural integrity monitoring 
applications has long been sought. Developments over the last several years have gone far 
towards obtaining a better understanding of the AE process in actual structural components 
(Refs. 1-5). These include the development of a model consisting of a repeatable source of 
known type, a structure for which the propagation of transient elastic waves can be computed, 
sensors whose transfer characteristics are either known or can be determined experimentally 
and the development of new signal processing algorithms for analyzing the detected signals. 

The associated problem, sometimes called the "inverse source problem", has as its aim the 
characterization of a source in terms of its location in the structure, its type, specified in 
terms of force or moment tensor components, its spatial distribution, if the source is extended, 
and its temporal or dynamic characteristics. This paper gives an overview of the work recently 
completed in our laboratory on aspects of this topic. 


THEORY AND WAVEFORM PROCESSING 


The displacement signals u^ detected at a receiver location, x, in a structure from an 
arbitrary source f(x',t) iocated at x' having source volume V can be written as 


u k = 


E(x',t) * G. k (x/x',t) dV 


(1) 


where the term Gj[< represents the Green's function of the structure. Insight into this 
formulation can be obtained if one makes a Taylor series expansion of the Green's function 
about the source point, 2 ', located at, say 0, to give 


u (x,t) = f.(0,t) * G.. (x/0,t) dV + x. f.(0,t) * G.. .(x/0,t) dV' 

*' j\i i J K J y 1 J JXj 1 

0 0 


( 2 ) 
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Tne first term in this equation, representing the monopole contribution, can be rewritten in the 
form 


or 



. f-(0,t) dV> 

v o' 


* Gj k (x/0,t) 


uj c m) = S.(t) * G jk (x/0,t) 


(3) 


The second term in Eq. (2), representing the dipole contribution, can be rewritten in a similar 
fashion as 


or 



{ x'. f .(x',t) dV' } 

Jv 0 ‘ 


* G jk>i (s/0,,) 



W) 


The term, M.. , called the moment tensor of the source, can be used to represent any 
combination o^double forces (Ref. 3) to model the source. It can be shown that the moment 
tensor can be written in terms of the crack's displacement jump and its area in a material (Ref. 
6). If each of the components of the moment tensor of the source have the same time 
dependence, then the spatial and temporal characteristics of the source are separable. Thus, 
the moment tensor can be written as a product of the moment density tensor, m.., and the and 
the source-time function, S(t), J 


M.. = S(t) *m.. (5) 

i) 

In recovering the characteristics of a dipolar source from the measured signals, u z , one uses 
knowledge of the Green's function G and inverts Eq. (4) to recover S(t) and the components 
m.. . In the most complicated deformation processes and for an arbitrary crack orientation, all 
six 1J components of the moment density tensor may need to be evaluated. As of now, this has 
not yet been done. Here we show the results obtained for several simpler cases. 

The simplest case is that of a point source of known type for which only the temporal 
aspects of the source are to be recovered from the detected signals. In that case, Eq. (3) can 
be rewritten as a simple convolution equation in time 

u(t) = S(t) * G(t) (6) 

This equation has been modified to include the transfer characteristics of a non-displacement 
transducer, whose response is band-limited. Recovery of the source function, S(t) , then 
involves two measurements, one of which is a calibration experiment using a known step 
excitation source applied at the source position. Examples of processing near -field and far- 
field signals from an axisymmetric source have been published previously (Refs. 7-9). 

For the case described by Eq. (4) in which both the time-function and moment-density 
tensor components of a source are to be recovered, several processing algorithms have been 
developed. One of these is based on an algebraic scheme in which the amplitudes of particular 
ray arrivals in the detected signals are directly related to the amplitudes of the Green's 
function (Ref. 10). Once these are determined, the source time-function can be recovered as 
before. The problems related to this scheme result principally from the amplitude 
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measurement precision required of specific features on a waveform. An alternative is a double 
iterative scheme demonstrated by Michaels and Pao (Ref. 9). The procedure is initiated with a 
guess of the moment tensor magnitudes from which a solution of the source time-function is 
obtained. This solution is used, in turn, to evaluate a new set of moment tensor amplitudes, and 
so on. The iteration process is continued until the results obtained in the iterations converge to 
a final answer. The procedure appears to be initial guess dependent and it may also converge 
slowly for certain configurations of the source moment tensor. It has been found that an 
appropriate initial guess for the double-iterative scheme can be obtained by making use of the 
results obtained from the algebraic processing scheme. 


RESULTS 

Demonstration of any signal analysis procedure for AE source characterization must begin 
with a controlled deformation process in a material. Reproducible, deformation-related AE 
sources whose location is known precisely, include indentation-produced cracks in brittle 
materials and fatigue-generated cracks in thin specimens and these were investigated in this 
work. 

a. Fatigue Cracks 

In detailed fatigue experiments, conventional test parameters were recorded in 
combination with electrical crack length measurements and a total of five AE waveforms, 
detected at points around the crack region. Two sample signals recorded equidistant from the 
crack front are shown in Figure 1. The four equivalent forces or moment tensor components by 
which an arbitrary dipolar source in a thin specimen can be modelled are shown schematically in 
Figure 2. It is found that the dynamic Green's functions corresponding to the M I? and M 7 . 
source components exhibit a non-zero signal amplitude between the P- and S-wave arrivals; 
From observations as this, each corresponding to an event, it is concluded that only two signals 
are needed to determine without ambiguity the characteristics of the source. It was found that 
some of the events could be adequately described in terms of only two double forces, i.e., M,. 
and M 22 • Other events required also a contribution from the M, 2 and M 2 , components. Table 
1 shows the results obtained for two events one of each type. More than thirty-eight signals 
were processed in a typical test. In Figure 3 is the recovered source-time function obtained 
with the double-iterative processing scheme. It appears that the crack displacement profile had 
a risetime of approximately 400 nsec. The risetime may actually be shorter than this, since the 
recovered value corresponds to the upper bandwidth limit of the transducer /amplifier system. 

b. Indentation Cracks 

The pressing of an indentor onto the surface of a flat plate results in a stress field 
configuration which depends on the shape of the indenter and the mechanical properties of the 
plate specimen. If the specimen is a brittle material, cracks of specific shape form when and in 
those regions where the stress field exceeds the fracture strength of the material. The dynamic 
force released in forming the crack is emitted as an AE signal. A typical epicentral AE signal 
obtained with a spherical-tipped indenter is shown in Figure 4(a). Typical Hertzian cone cracks 
are obtained when such indenters are pressed against a specimen. It has been shown that the 
associated dynamic source is equivalent to a vertical axisymmetric force acting on the 
structure (Ref. 11). On unloading, one or more lateral vent cracks may form whose associated 
AE signals exhibit completely different characteristics as shown in Figure 4(b). In some cases, 
these may be adequately modelled by a single horizontal dipole acting near the surface of the 
plate. This appears to be the case, when only one lateral vent crack forms. The detected AE 
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signals resemble those obtained when a conical indenter is pressed against the specimen. In the 
latter situation, a single penny-shaped crack is formed. Evaluation of the stress field under the 
conical indenter, obtained from the Boussinesq solution (Ref. 12), shows regions of tensile and 
compressive hoop stresses. The crack is generated in the region of tensile hoop stresses whose 
contours resemble a penny. 

Experiments were carried out with an eight-channel AE system consisting of point 
piezoelectric transducers with a 1.2 mm active area situated at various angles at a distance 
approximately 2h (h - plate thickness) from the indenter. Sensors whose response was nearly 
uniform in amplitude were selected for the array. The signals were amplified by 60 dB with 
preamplifiers whose bandwidth ranged from 10 kHz to 2 MHz. They were subsequently 
digitized in transient waveform recorders operating at a sampling rate of 30 MHz with 10-bit 
resolution. 

Evaluation of the moment tensor components of the penny-shaped crack via one of the 
procedures, listed in the previous section often did not yield reliable results, as they were 
strongly initial guess dependent. It was clear, however, when comparing the detected 
waveforms to synthesized AE signals, that the model of a tensile crack given in terms of three 
orthogonal dipoles of strength, X + 2 y , X and X , where X and u are the Lame constants, 
may not be appropriate to model an indentation-produced crack in a small specimen as 
described here. 

An alternative approach for recovering the ratios of the moment tensor components 
involves measurement of the radiated field of the AE signal from the source. The radiated 
amplitude measured at a distance R from a dipolar source, can be written in terms of the 
displacement amplitude, u z , and this is expected to be of the form 

u z (R, t) = aj(t) + a 2 (t) cos 2 0 (7) 

where 0 is the orientation of the normal to the crack. For a single horizontal dipole source, 
aligned parallel to the crack normal, the term a , (t) is expected to be near zero. For three- 
dipole crack model, its value is expected to be approximately 0.28. Eq. (7) can be rewritten to 
account for an arbitrary crack orientation relative to the sensor array 

u z (R, t) = a^(t) + a 2 (t) cos 2 ( 0 - 0 C ) (8) 

where 0 is the orientation of the normal of the crack relative to the orientation of the sensor 
array. c 

By a proper choice of angular distribution of sensors around the indenter, an eight-channel 
AE system is adequate to determine the field characteristics of the indentation-generated 
cracks. Each channel of the entire system was calibrated simultaneously by breaking a glass 
capillary 0.05 mm in diameter at the crack position on the specimen. The result was an 
axisymmetric step excitation source. The amplitudes of both the first P- and S-wave arrivals 
were measured. The corresponding amplitudes in the crack-generated signals were normalized 
with the calibration amplitudes. The result obtained after one indentation is shown in Figure 5. 
The solid and dashed lines correspond to the results produced from a non-linear least-square fit 
of Eq. (8) to the measured data. From the curve fit, the parameters of the radiated field can be 
extracted. Once these are .known, the measured amplitudes at each transducer can be reflected 
across the various symmetry planes to generate image data points with the results obtained 
shown in Figure 6. Inspection of the amplitude ratio a./a 2 shows that the tensile cracks 
produced by the conical indenter are best modelled as 1 resulting principally from a single 
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horizontal dipole. The orientation of the crack, relative to the array configuration, is 
recovered immediately when Eq. (8) is fit to the amplitude data. The difference between the 
observed and recovered crack orientation angle from seven indentations was always less than 
6° and on average less than 3°. 

The radiated field of the signals emitted from lateral vent unloading cracks was also 
measured. It was found that the formation of multiple cracks results in a radiated field which 
radiate more uniformly in all directions and clearly, the radiation does not resemble that of a 
single crack. 

It is possible to process an epicentral displacement signal with the double-iterative 
processing scheme to recover the individual magnitudes of the moment tensor as well as their 
time function. Figure 7 shows the recovered source time function for an indentation crack. 


EXTENDED SOURCE 

Both the forward and inverse problems of the extended source have been solved using 
synthetic data. A series of examples dealing with two source types and several linear 
distributions on the surface or through the interior of a plate specimen have been computed 
(Ref. 13). In Figure 8(a) is shown the case of a distribution of horizontal dipole sources, 
extending to 0.75h into the interior of the plate and whose time function was a Heaviside step 
function. It has been proposed that such a source may be used to model the interaction of 
intense synchrotron-generated pulses of x-rays striking a specimen and penetrating into its 
interior (Refs. 14-15). Figure 8(b) shows a comparison between the signal from a single dipole 
very close to the top surface of the specimen and the signal synthesized from forty source 
points whose strength distribution along the z-axis was given by e-0.5z. After only ten 
iterations in the double-iterative processing scheme, the input source function was recovered 
from the synthetic input data within a percent and the corresponding recovered source time 
function was indistinguishable from the input. 


CONCLUSIONS 

Acoustic emission signals detected from cracks formed during indentation and fatigue 
loading experiments have been processed to recover the moment tensor components of the 
crack source and their time function. Several signal processing schemes were demonstrated 
including one based on ray amplitude measurements, a double-iterative processing scheme and 
an AE radiation pattern analysis method. A signal processing algorithm was also demonstrated 
with synthetic AE waveforms corresponding to an extended source of emission by which the 
signals can be processed to recover both the spatial extent of the source and its time function. 
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CHARACTERIZATION OF FATIGUE CRACKS IN 7075 -T6 


Uj, = CnGriU +C22 .Gt22.,2. + C 12 . (kx,z + C 21 

Event #3^- ; Model : Two Doable Forces 
Al^ebrO-tC J C 11 I = | : 0,258 

Doable Iterate. on ? Cn ' C 22 = I 1 0.266 

Event #35 Model "Two Doable Forces, Two Single Coaples 
Algebraic ; C,i C 22 : C 3 j ‘.Ch* 1 : 0.324 : -0.271 : 0.237 
Doable I ter at Con ; Cu : C 22 :Gi : C*=W9€ :o3oz :-oa7\'oao3 


Table I 



VOLTAGE *10 


FATIGUE CRACK SIGNALS IN 7075-T6 


FATIGUE CRACK: SENSOR #1 FATIGUE CRACK: SENSOR *2 



TIME (useo) TIME (used 

Figure 1 


EQUIVALENT FORCES / MOMENT TENSOR COMPONENTS 



Figure 2 
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RMPLITUDE *1Q 


RECOVERED SOURCE-TIME FUNCTION 



Figure 3 


UNLOADING CRACK: I4.5NT/MAX: 35.3NT 


120° CONICAL INDENTER - DIPOLE RADIATION 


90 



P - Wave 
S - Wave 


270 

90 


Figure 5 









RECOVERED SOURCE -TIME FUNCTION (CRACK 9) 



Time (^tsec) 

Figure 7 


DIPOLE SOURCES: EXPONENTIAL DISTRIB 



NORMALIZED TIME t = CpU/h 


Extended dipole source model Synthetic normal displacement signals 

corresponding to several source 
distributions 

Figure 8(a) Figure 8(b) 
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